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@ Procede de reconstruction d'images tridimensionnelles d'un objet par des mesures utilisant un 
rayonnement conique et un reseau bidimensionnel de detecteurs. 

@ Proced§ de reconstruction d'images tridimensionnelles d'un objet (3) dans leque) la reconstruction 
est obtenue en calculant des sommes de projections sur des plans extraits d'un rayonnement 
d'attenuation conique traversant I'objet ou d'une emission provenant de I'objet suivant une collimation 
conique. Le reseau bidimensionnel de detecteurs (4) utilise a cet effet est a une inclinaison (Q 
quelconque par rapport a Taxe de rotation definissant le repere dans lequel I'image est reconstruite, par 
suite d'un choix volontaire. On donne les formules qui permettent d'obtenir une reconstruction a faible 
distorsion apres avoir mesure les parametres de deplacement du reseau et du faisceau de rayonnement 
conique. Les erreurs de position qui resultent de defauts de construction sont 6galement prises en 
compte dans les calculs apres avoir 6te mesur6es. Un protocole de mesure des parametres et ie 
procede de calcul associe sont decrits. 
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L'invention se rapporte a un precede de reconstruction d'images tridimensionnelles d'un objet a I'aide de 
mesures utilisant un rayonnement conique et un reseau bidimensionnel de detecteurs. 

Elle peut etre consid6r6e comme un perfection nement des inventions decrites dans les brevets precedents 
attribues a la demanderesse (brevet europeen 0 29 2402 et brevets francais de numeros d'enregistrement 90 
14957 et 90 14958). 

Un precede de reconstruction de I'image d'un objet tridimensionnel a partir d'une serie d'irradiations par 
un faisceau conique tournant autour de I'objet et qui vient frapper un ecran couvert d'un reseau bidimensionnel 
de detecteurs est plus precisement decrit dans le premier de ces brevets : un algorithme utilisant Inversion 
de ce qu'on appelle la derivee premiere de la transformee de Radon de i'attenuation du rayonnement du fais- 
ceau a travers I'objet a reconstituer est utilised 

Les calculs exposes la ne sont cependant plus valables quand le reseau de detecteurs est maJ place par 
rapport a ia trajectoire, ou quand le point focal est d6plac6 de sa position theorique, car la ligne de projection 
orthogonale du pointfocal surle reseau de detecteurs (I'axe focal) ne passe plus par I'origine du repere d'objet 
Une grande partie des calculs est en effet faite dans ce repere, qui est choisi sur I'axe de rotation du point 
focal et du reseau de detecteurs quand ceux-ci parcourent des trajectoires circulaires et dans le plan des tra- 
jectoires quand celles-ci sont coplanaires. Or ces situations sont extremement frequentes car elles peuvent 
provenir de defauts d'alignement et de position nement des mecanismes porteurs de la source et de recran, 
defauts qu'on n'a pas toujours le loisir de pouvoir corriger. 

Parfois ces situations sont meme volontairement recherch6es : un exemple classique est celui de I'examen 
de I'interieur du crane ou le reseau de detecteurs parcourt une trajectoire situee sur un cercle dont I'axe est 
parallele a I'axe du patient. 

On desire utiliser I'appareil pour examiner des regions internes a des hauteurs differentes. Pour des ran 
sons de cout et de simplicite, il est avantageux de pourvoir le reseau de detecteurs d'un collimateur unique 
pour toutes les situations d'examen. 

Or, les distorsions de la fonction de transfert de modulation qui accompagne les mesures et la reconstruc- 
tion, et done les distorsions de I'image, sont moins grandes dans le plan de la trajectoire du point focal, ce qui 
explique qu'on cherche a faire couper la region a examiner par ce plan. 

On est alors amene a vouloir regler rinclinaison du reseau de detecteurs sur le crSne du patient, en le dis- 
posal obliquement de facon que le point focal suive une trajectoire circulate dont la hauteur est reglee en 
consequence. De plus les detecteurs suivent ainsi au niveau des epaules des trajectoires de plus grand dia- 
metre que ceux qui sont au niveau du sommet du crane, ce qui permetd'enfoncer la tete plus profondement 
dans le dispositrf et ainsi d'abaisser le plan de la trajectoire du point focal en direction du cou. Si le point focal 
correspond a une source materieile de rayonnement, il suff it de la relier au reseau de detecteurs par une 
connexion mecanique rigide. 

Les calculs de reconstruction des images sont sensiblement plus compliques dans ces conditions geo- 
m&riques plus generales. Certaines publications les mentionnent ou les evoquent, au moins dans certains cas 
particuliers. II est ainsi decrit comment proceder quand la reconstruction est effectuee par ralgorithme de Fel- 
dkamp dans lequel le faisceau conique est decompose en un empilement de faisceaux plans en forme d 'even- 
tail (article de Manglos, Jaszczak et Greer "Cone-beam SPECT reconstruction with camera tilt", paru en 1 989 
dans "Physics in Medicine and Biology", n°34(5), p.625 et abrege de Cao et Tsui paru en 1 991 dans le "Journal 
of Nuclear Medicine", n°32(5), p.1 066) ou par un algorithme de calcul iteratif dans lequel ia contribution de che- 
que point a I'attenuation du rayonnement est calculee a partir d'une valeur initiate en calculant a chaque etape 
des facteurs de correction deduits du rapport entre les mesures et les projections de I'estimation de ia fonction 
obtenue a I'etape precedente (article de Manglos, Jaszczak et Mac Afee "Maximum likelihood reconstruction 
for cone-beam SPECT with camera tilf paru en 1 989 dans "IEEE Transactions on Nuclear Science", n°36(1). 
p. 11 17), mais rien ne semble avoir ete effectue dans le cas ou I'algorithme implique I'inversion de la transfor- 
mee de Radon de la fonction ou de sa deYivee premiere. 

L'invention comble cette lacune, aussi bien dans le cas evoque jusqu'a present d'un faisceau conique emis 
une source ponctuelle externe et qui traverse I'objet en etant attenue que dans le cas parfaitement equivalent 
d'un reseau bidimensionnel de detecteurs collimates vers un pointfocal unique derriere I'objet et sensible au 
rayonnement d'un produit emetteur de rayonnement qui impregne I'objet, comme un marqueur contenu dans 
un liquide injecte au patient et qui se r6pand dans un organe a explorer apres avoir ete transports par le reseau 
veineux. II faut aussi observer que l'invention peut s'appliquer a un nombre quelconque de rSseaux de detec- 
teurs qui contribuent chacun a une partie de mesure et sont deplaces de quantites differentes. 

L'invention concerne done sous sa forme la plus generate un proc6d6 de reconstruction d'images tridi- 
mensionnelles d'un objet defini par des valeurs prises par une fonction sur des points de I'objet, la fonction 
etant une propriete d'un rayonnement conique ayant un point focal et passant a travers I'objet, dans lequel la 
fonction est calculee par I'intermediaire de sommes de la fonction sur des plans parametres passant par au 
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moins un point de I'objet et definis dans un rep£re d'objet comportant une origine, le proc6d6 comportant au 
moins une s6rie de mesures prises dans differentes positions du point focal autour de I'objet, chaque somme 
de la fonction sur lesdits plans parametres dans le repere d'objet etant calculee a partir d'au moins une des 

5 mesures, chacune des mesures etant men6e avec un reseau bidimensionnel de detecteurs du rayonnement 
orientes vers le pointfocal, etqui estd6plac6 sur une trajectoire circulaire autour de I'objet perpendiculairement 
a un axe de rotation, le point focal etant d£plac6 en m6me temps que le reseau de detecteurs de maniere a 
se projeter orthogonalement sur un point fixe du r6seau des detecteurs, un algorithme d'inversion etant utilise 
pour retrouver la valeur prise par la fonction sur les points de I'objet a I'aide des sommes de la fonction sur les 

10 plans parametres, I'algorithme comprenant des etapes de determination d'au moins une position du point focal 
appartenant a chacun desdits plans parametres, caracteris6 en ce que I'axe focal, d6fini comme ta ligne de 
projection orthogonale du point focal sur le reseau, passe a l'6cart de I'origine et en ce que ies etapes de d6- 
termination comprennent des calculs ou interviennent des parametres exprimant l*6cart entre I'axe focal et I'ori- 
gine et I'inclinaison entre I'axe focal et I'axe de rotation, cette inclinaison pouvant £tre quelconque. 

15 La geometrie d'acquisition est d6crite d'une part par un jeu de parametres intrinseques a savoir la distance 

focale FGd, la position Gd de la projection orthogonale du point focal F sur le reseau de detecteurs, les pas 
d'echantillonnage s6parant les detecteurs sur le reseau, et d'autre part par un jeu de parametres extrinseques 
servant a positionner le point focal et le reseau de detecteurs consideres par rapport a I'axe de rotation, a savoir 
la position du point Od projection conique de I'origine 0 sur le reseau de detecteurs, la distance FG s6parant 

20 le pointfocal et le plan parallele au detecteur et contenant I'origine 0 du repere d'objet, I'inclinaison du vecteur 
perpendiculaire aux detecteurs par rapport a I'axe de rotation, I'angle de rotation des lignes du r6seau de de- 
tecteur par rapport a la droite d'intersection du plan de rotation du point focal avec le detecteur, la position an- 
gulaire de la premiere acquisition. 

On va maintenant decrire i'invention plus en detail a I'aide des figures suivantes annex6es a titre illustratif 

25 et non limitatif : 

- la figure 1 repr6sente un exemple d'utilisation de I'invention ; 

- et les figures 2 a 7 represented les elements g6ometriques mis en oeuvre et commentes dans la suite 
du texte. 

La figure 1 repr6sente done un patient 1 etendu sur une table 2 et dont la tete 3 est examinee a I'aide d'une 
30 camera 4 de rayons gamma, qui tourne autour de la tete 3 selon un axe de rotation 5 etabli dans le sens de 
la hauteur du patient. La camera 4 estformee essentiellement d'un r6seau plan (ou plus g6neralement bidi : 
mensionnel) de detecteurs coliimates de maniere que les rayon nements qu'iis recoivent convergent tous vers 
un pointfocal F. Elle est fix6e a une structure m6canique circulaire non representee - mais dont les brevets 
cites fournissent des exemples - de facon que le point focal F est solidaire de la camera 4 et parcourt une tra- 
35 jectoire T en forme d'orbite circulaire perpendiculaire a I'axe de rotation 5 et dont le centre 0' appartient a cet 

axe. La trajectoire T s'etend ici a hauteur du cerveau du patient 1 . La camera est inclin6e d'un angle ( | - U 

par rapport a I'axe de rotation 5 ; I'angle £ sera d6f ini plus pr6cisement d'ici peu. II est suppose 6gal a | dans 

40 les brevets anterieurs ; I'appareillage associe a cette invention est de preference constitue de facon a permettre 
de modifier Tangle £ a volonte pour d6placer la trajectoire T et la placer a hauteur de la region de la tete 3 pour 
laquelle on manifeste de IMnterfit, car une plus grande nettete est obtenue a hauteur de cette trajectoire T. 
Les points de I'image du patient sont d6f inis par leurs coordonn6es dans un repere d'objet qui est un repere 

^ cartesien d'origine 0 place sur I'axe de rotation 5 et def ini par trois vecteurs directeurs : le vecteur ? oriente 

vers la gauche du patient 1, le vecteur J oriente vers I'avant et le vecteur "jT oriente vers le bas et colin6aire 
a I'axe de rotation 5. Le but de I'invention est de perfection ner des algorithmes existants de reconstruction de 
I'image a partir des vues prises par la camera 4 pour les rendre utilisables dans ce cas g6n6ral. Les relations 

^ g6ometriques (angles et distances) entre les deux series d'eiements sont supposees connues par une mesure 
anterieure du type de celle pr6sent6e dans la suite du brevet, quelle que soit leur origine, d6faut de construction 
de la camera 4 ou de son m6canisme de support ou r6sultat d'un choix volontaire de I'inclinaison de I'axe focal 
du r6seau de detecteurs. On ne tiendra pas compte dans cet expose des phenom6nes parasites tels que I'auto- 
att6nuation du rayonnement a I'interieur de Porganisme examine. 

^ Nous nous placons dans le cadre des methodes analytiques. La cartographie de I'activite d 'attenuation 

ou d'emission est decrite par une fonction tridimensionnelle f representant un ensemble de coupes transverses 
paralieies. Pour chaque position du pointfocal F, chaque detecteur 6 est relie au pointfocal F par une droite 
de projection unique. La projection en geometrie conique Xf associe a chacune de ces droites une integrate 
de la fonction f . La reconstruction revient a r6soudre le probieme inverse, e'est-a-dire a mettre en oeuvre i'en- 
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chamement des transformations necessaires pour calculer f a partir de Xf. L'outil d'analyse est la transforma- 
tion de Radon 3D, def inie comme I'ensemble des valeurs integrales sur les plans de I'espace. Dans ia pratique, 
on se restraint aux plans qui rencontrent I'objeL On associe a chaque plan P son point caracteristique C, pro^ 

jection orthogonale de I'origine 0 sur ce plan, et les coordonnees spheriques ( p, n* ) de ce point C. 
Le point C, def ini par ses coordonnees x, y et z dans le repere d'objet 



peut egalementetredefini en coordonnees spheriques par la relation dc-pn ou "ft est un vecteur unitaire 
et p un rayon. II faut alors introduire les coordonnees telles que la longitude 9 et la colatitude 9 pour def Inir 
I'orientation du vecteur "ft . 

15 Pour une position donnee du point focal F et pour un plan passant par ce point, nous avons etabli une for- 

mule exacte, reliant la projection en geometrie conique Xf, a la derivee premiere de la transformed de Radon 
R'f. L'operation principale consiste a sommer le long de ia droite d'intersection du plan avec ie d&ecteur les 
derivees premieres de la projection. Ceci permet d'affecter une valeur a tous les plans passant par cette po- 
sition du plan focal. Apres avoir procede ainsi pour chaque incidence de mesure, I'operation de rearrangement 

20 consiste a redistribuer ces valeurs integrales associees a chaque plan, du systeme de coordonnees lie aux 
positions du point focal, sur le systeme de coordonnees spheriques du domaine de Radon. Pour retrouver la 
fonction f, il reste a inverser la derivee premiere de la transformee de Radon comme il est d6crit dans le brevet 
mentionne en premier. L'algorithme procede en deux etapes, enchainant d'abord des retroprojections bidimen- 
sionnelles en geometrie parallele pour chaque meridien du domaine de Radon, puis des retroprojections bidi- 

25 mensionnelles en geometrie parallele suivant les plans axiaux de I'objet. 

La trajectoire du point focal est usuellement circulaire. En caracterisant Information acquise par les plans 
rencontrant la trajectoire, notre analyse met en evidence l f existence d'une zone d'ombre dans le domaine de 
Radon, liee aux plans qui rencontrent I'objet mais pas la trajectoire. Cette zone d'ombre est une source d'ar- 
tefacts sur I'image reconstruite. Ceci explique pourquoi il est souhaitable que le plan focal soit centre sur ia 

30 region a reconstruire. Notre algorithme permet d'attenuer ces distorsions en comblant cette zone d'ombre par 
interpolation. 

Pour simplif ier les calculs, on effectue un parametrage sur un plan de detection Pdet, parallele au plan 
des detecteurs Pm auquel les detecteurs 6 appartiennent, et qui passe par I'origine 0. Les calculs menes sur 
le plan de detection Pdet sont valables quelle que soit la distance du plan des detecteurs Pm, a un facteur de 
35 grossissement pres. On designe par A un point quelconque du plan de detection Pdet, par G la projection otho- 
gonale du point focal F sur le plan de detection Pdet, et par Od, Ad, Gd les points du plan des detecteurs Pm 
tels que F,0 et Od, d'une part, F,A et Ad ensuite, et F,G et Gd d'autre part soient alignes ; mais dans le cas 
present, G et O sont distincts contrairement aux brevets precedents. Les calculs menes sur le plan des detec- 
teurs Pm ou, de maniere equivalents sur le plan de detection Pdet utilisent un repere d'origine 0 dont deux 
vecteurs ^ et ^ sont disposes paralielement aux deux directions principales (lignes et colonnes) du reseau 

de detecteur 6 et dont le troisieme "3 est dirige vers le point focal F dans une direction perpendiculaire Ce 
repere est nomme le repere d'acquisition. 

Trois angles de rotation peuvent etre definis (figure 4) pour decrire la position respective des reperes d'ob- 
jet et d'acquisition. L'angle de precession y autour de I'axe ol transforme le repere 



en un repere intermediate 

II est associe a la rotation de la camera 4 et du point focal F sur la trajectoire T. 
L'angle de nutation ^ fait passer du repere 

<fu-?'t 

au repere 
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et est done une rotation autour de I'axe def ini par le vecteur if ' ; il exprime TincNnaison de la camera 4 par 
5 rapport a I'axe de rotation 5. 

L'angle de rotation propre 8 est une rotation autour de I'axe def ini par le vecteur y et I'origine 0 et qui fait 
passer du repere 

10 o u 1 v 1 w 

au repere 

Ouvv? ; 

15 

il exprime I'inclinaison des itgnes des detecteurs 6 de la camera 4 par rapport a ('intersection du plan de la 
trajectoire T avec le detecteur. 

U ne connaissance precise de 8 et surtout de £ est necessaire pour eviter les imprecisions de calcul pendant 
la reconstruction. La connaissance exacte de y est utile pour connaTtre avec precision la position exacte de 
20 I'image par rapport a I'objet, mais les imprecisions sur cet angle n'ont pour consequence qu'une erreur d'orien- 

tation de I'image reconstruite. Les vecteurs ^ , "v* et "w . Q ui sont li6s £ ,a camera 4 etaux plans des detecteurs 

Pm et de detection Pdet, sont illustres aux figures 5 et 6 : u* et "v appartiennent a ces plans et sont orientes 

25 respectivement en direction des lignes et des colonnes du detecteur 6, et ^ est orthogonal a ce plan et dirige 
vers le point focal F. Le point O est I'origine des coordonnees de ce repere. 
L'expression de la transformee de Radon est definie par la formule (1) : 



30 



Rf 



MeP(p.n ) 



f(M) dM 



(1) 



celle de sa derivee par la formule (2) : 



R*(p.Tf)-«*£ '( p.Tf) ,„ 



40 

et les expressions de leurs inversions par les formules (3) et (4) ; 



50 



55 

Lafonction qui fait i'objet de la transformee est designee part la transformee elle-meme par Rf, M designe 
les points du repere d'objet et S 2 designe I'ensemble des vecteurs unitaires de I'espace. 

Si maintenant nous appelons A un point quelconque de Intersection d'un plan P (p, "ft ) et du plan de 
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detection Pd6t, et si nous appelons et n les vecteurs deduits des vecteurs u et parune rotation 

d'angle a autour du vecteur J pour que le vecteur t " soit directeur de la ligne qui constitue cette intersec- 
tion, la formule (5) : 



10 



15 



RT OF "I? *Tv I l,GF,, f 3Yf , 

Rf ,Or. n.n j = _ J g^r (S 1 A(q»)) dq" (5 ) 

II G FA n* II H 

exprime la valeur de la derivee premiere de la transformee de Radon, p" et q" etant les coordonnees du point 
A dans le repere O 

u"v M 



(OA-p"u w +q"?") et Yf designe la projection Xf de la fonction mesuree sur le rayon qui mene du point focal F 
20 au point A apres une pond6ration telle que 



25 



Yf (F, A) = IlfjlL xf(F A) 

\\fk\\ ' A) 



!ci, O et G sont distincts, le vecteur o*p est ega! a 



30 



pgu+qg^v+FGv, 



35 



40 



soit(x'F) u "+(y'F) j '+(zF)k. Les definitions precedentes nous permettent d'obtenir le systeme d'equation (6) : 



y f = cosi . sin5 . po +• cosq . cos5 . qc -sinq . FG 



(6) 



= sine . sin5 . pc, + si nq . cos8 . qc + cos^ . FG 



45 et si A\j/ est Tangle entre et 6*F \ ou F est la projection orthogonale du point focal F sur le plan def ini par 



I'origine 0 et les vecteurs if et ?, on peut dSduire assez facilement les equations (7) et (8) pour une trajectoire 
circulaire du point focal F, d'autres formules pouvant etre deduites pour d'autres formes de la trajectoire T : 



50 



V|/ = ip - Ay +■ arc sin 



p-zr : eos8 
OF . sine 



(7) 



55 



V = <p + it - Av • arc sin 



-z F ■ cos9 l 
)F' . sine J 



P 

OF 



(8) 
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En remarquant que C( p, 9, q>) est identique a C( p, -G, q> +n), on peut se restreindre a n'utiliser que Tune 
d'entre ces equations, par exemple la (7). 

Ces equations (7) et (8) permettent d'obtenir pour chaque point du volume caracteristique de I'objet d6f ini 
par ses coordonnees spheriques p, 0 et <p !es valeurs des angles de precession vy pour lesquels les mesures 
contenues dans le plan associe, et done la derives premiere de la transformee de Radon en ce point peut etre 
obtenue. Ces angles sont generalement au nombre de deux puisque les plans sont generalement secants en 
deux points a une trajectoire circulate. 

Les interpolations pour deduire les positions effectives du point focal F qui vont servir au calcul - il s'agit 
des positions voisines de cette position obtenue par calcul en utiiisant ies equations (7) et (8) et pour lesquelies 
des mesures ont ete faites - ainsi que les r^tro- projections pour reconstruire ('image en des points exprimes 
dans le repere d'objet a partir des sommes calculees sont identiques a celles du premier brevet 

II faut toutefois examiner comment il convient d'effectuer les interpolations pour combler la "zone d'ombre" 
des mesures avant d'accompiir les calculs d'inversion. On se reporte pour cela a la figure 7 : il est possible 
de demontrer facilement que, si la trajectoire du point focal F est un cercle, les points caracteristiques C as- 
socies a des plans sur lesquels la somme de la fonction f peut etre mesuree sont inclus dans un tore To obtenu 
par la revolution d'un cercle de diametre O'F autour du centre O' de la trajectoire T. Si on designe par p\ 9 et 

<p les coordonnee s spheriques du point C a partir du centre 0\ avec o"^c = p' "rt v et s ' 6"V = zF t , II est 
possible d'exprimer la relation entre p' et par la formule : 

p' = p - zF. cosG. (9) 

Le principe adopte dans le brevet mentionn6 en premier (la difference par rapport au premier brevet est 
que 0 et 0' sont distincts) (interpolation d'ordre zero) consiste a associer a un point C de la zone d'ombre la 
valeur mesuree sur un point d'une paire (O ou C+) a la limite de la zone d'ombre, de meme rayon p' et de 
meme longitude q> ; les colatitudes de ces deux points sont +0i et -9i. On choisit le point dont la colatitude est 
la plus proche de 0, ici C+. 

II suff it alors d'assimiler ia valeur R'f (p\ 9, <p ) du point C def ini par ses coordonnees p', 9 et <p a partir de 
I'origine 0' a la valeur mesuree R'f ( pi, Gi, pi) ou : 
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pi = P % + zF 



1 - 



OF' 2 ) 



1 /2 



(10) 
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9i = ± arc sin 



(---) 

40 \ OF' J 

selon que C+ ou C- est choisi, 
et 

45 P' = P 

Tous ces calculs ne sont cependant accomplis avec succes que si les parametres geometriques sont 
connus a Tavance avec une precision suff isante pour decrire correctement le systeme d'acquisition des me- 
sures. 

Certains parametres, a savoir les coordonnees du point Gd sur le plan de detecteurs Pm, la distance focale 
50 FGd et les distances d'echantillonnage separant les points d'acquisition sur le detecteur peuvent etre evalues 
en disposant une grille ou un dessin analogue devant le plan des detecteurs, et a I'illuminer par une source 
plane qui peut etre un petit reservoir rempli de liquide radioactif af in de recueillir I'image de ia grille sur le plan 
des detecteurs. La grille est integree a un mecanisme qui permet de la placer parallerement au plan des d6- 
tecteurs et a plusieurs distances. L'image se dilate et se contracte done sur le plan des detecteurs Pm suivant 
55 qu'on en approche ou qu'on en eloigne la grille. Le point Gd est le centre des dilatations et des contractions, 
le grossissement peut etre evalue en fonction des distances de ia grille au plan des detecteurs Pm et a la source 
et de la taille de l'image de ia grille. 

Un dispositif de ce genre est decrit dans I'article de Chang, intitule "New Methods of Examining Gamma 
Camera Collimators ,, et paru en 1988 dans "The Journal of Nuclear Medicine", vol.29, n°5, p.676 a 683. 
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On peut remplacer la source lumineuse unique et la grille par un reseau de sources aux noeuds d'une grille 
virtuelle. La plaque qui porte les sources est placee a des distances variables des detecteurs 6 par un jeu d'en- 
tretoises. 

Les angles <j>, \ et 8 et la distance FG peuvent etre obtenus par une serie d'experiences, en partie decrites 
dans un article publie ("Estimation of geometrical parameters and collimator evaluation for con-beam tomo- 
graphy", par Gullbera, Tsui, Crawford, Ballard et Hagius, paru en iggo dans Medical Physics n°1712) selon 
lesquelles on dispose au moins une source de rayonnement fixe et ponctuelle et on fait tourner le reseau de 
detecteurs pour obtenir une serie d'images de la source. La comparison de ces images en liaison avec les 
endro.ts d'ou elles ont ete prises permet de determiner les trois angles et d'en deduire certaines indications 
sur la quahte de collimation. Un algorithme de minimisation des incertitudes ou des erreurs selon une methode 
des moindres carres peut etre utilise. 

Enf in. si on utilise plusieurs detecteurs, I'ecart entre les detecteurs est donne par la difference entre les angles 
de precession y evalues selon la methode precedente pour une meme position des sources ponctuelles. 

L'invention peut etre appliquee tres commodement a des cas ou le point focal F et la camera 4 sont mus 
sur des trajectoires circulaires, mais elle peut aussi etre appliquee a des systemes a plusieurs trajectoires 
d acquisition comme dans les deux brevets francais cites. Enf in, on peut I'utiliser non seulement dans les al- 
gorithms utilisant la transformee de Radon de la fonction qui sert a la reconstruction de I'image ou la derivee 
prem.ere de cette transformee, mais encore dans les algorithmes qui component une transformee de Hilbert 
de la transformee de Radon ou de sa derivee premiere. 



Revendications 

1. Precede de reconstruction d'images tridimensionnelles d'un objet (3) def ini par des valeure prises par une 
fonction (f) sur des points (M) de I'objet, la fonction etant une propriete d'un rayonnement conique ayant 
un point focal (F) et passant a travers I'objet (3), dans lequel la fonction est calculee par I'intermediaire 
de sommes de la fonction sur des plans (P) parametres passant par au moins un point de I'objet et def in is 
dans un repere d'objet (O, ?,?.*) comportant une origine (0), le precede comportant au moins une 
serie de mesures prises dans differentes positions du point focal autour de I'objet, chaque somme de la 
fonction sur lesdits plans parametres dans le repere d'objet etant calculee a partir d'au moins une des 
mesures, chacune des mesures etant menee avec un reseau bidimensionnel (4) de detecteurs (6) du 
rayonnement orientes vers le point focal, et qui est deplace sur une trajectoire autour de I'objet perpen- 
diculairement a un axe de rotation, le pointfocal etant deplace en meme temps que le reseau de detecteurs 
de maniere & se projeter orthogonalement sur un point fixe (Gd) du reseau des detecteurs, un algorithme 
d'mversion etant utilise pour retrouver la valeur prise par la fonction sur les points de I'objet a I'aide des 
sommes de la fonction sur les plans parametres, I'algorithme comprenant des etapes de determination 
d'au moms une position du point focal appartenant a chacun desdits plans parametres associes, carac- 
terise en ce que I'axe focal, defini comme la ligne de projection orthogonale (FGd) du point focal sur le 
reseau passe a l'6cart de I'origine (O), et en ce que les etapes de determination comprennent des calculs 
ou mterviennent des parametres (A, <p, % ) exprimant l'6cart entre I'axe focal et I'origine et I'inclinaison entre 
I'axe focal et I'axe de rotation. 

2. Precede suivant la revendication 1, et pour une trajectoire plane (T) du point focal (F), caracterise en ce 
que les parametres exprimant l'6cart entre I'axe focal et I'origine comprennent la position de la projection 
orthogonale du pointfocal (F) surun plan de detection paralieie au reseau de detecteurs et passant par 
I'origine (O) du repere d'objet 

3. Precede suivant la revendication 2, caract6ris6 en ce que les calculs des etapes de determination des 
plans parametres comportent au moins I'une des formules (7) et (8), la trajectoire du point focal etant cir- 
culaire. 

4. Precede suivant la revendication 1 a 3, caracterise en ce que les etapes de determination comprennent 
des calculs ou interviennent ('angle de rotation (8) de lignes des detecteurs du reseau bidimensionnel par 
rapport au plan 



parallele a la trajectoire. 
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Proc6de suivant i'une quelconque des revendications 1 a 4, caract6ris6 en ce que I'algorithme calcule 
pour des plans passant par I'objet et rencontrant la trajectoire (T) du point focai en au moins un point ie 
taux de variation de la somme de la fonction sur ce plan a partir des mesures realises en ce point. 

Precede" suivant Tune quelconque des revendications 1 a 5, caracteYise en ce que ie resultat des sommes 
de fonction comprend une transformee de HilberL 

Proc£de suivant I'une quelconque des revendications 1 a 6 dans lequel Ie reseau plan de d6tecteurs tour- 
ne autour de la tSte d'un patient (1), caracte>is6 en ce que Ie r6seau des detecteurs est inc!in6 vers Ie 
corps du patient et des r6glages m6caniques sont accompiis pour que Ie point focal (F) suive une trajec- 
toire passant par une region d'interSt autour de la tste du patient. 

Proc6d6 suivant I'une quelconque des revendications 1 a 7, caracteris6 en ce que Pinclinaison entre I'axe 
focal et I'axe de rotation est quelconque. 

Proc6d6 suivant I'une quelconque des revendications 1 a 8, caract6ris6 en ce qu'il comprend une 6tape 
de complement d'un ensemble de valeurs (R'f) associ6es a des points de I'objet et obtenues par des cal- 
culs effectues sur les sommes de la fonction surles plans param6tr6s par des valeurs associees a d'autres 
points (C) de I'objet et obtenues en interpolant des valeurs obtenues par les calculs et associees a des 
points (C+, C-) de I'objet choisis en fonction des coordonn6es desdlts autres points par rapport a une in- 
tersection (O') entre un plan dans lequel la trajectoire du point focal est comprise et I'axe de rotation. 
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